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Abstract We discuss the degree to which radio propagation measurements diag- 
nose conditions in the ionized gas of the interstellar medium (ISM) . The "signal 
generators" of the radio waves of interest are extragalactic radio sources (quasars 
and radio galaxies), as well as Galactic sources, primarily pulsars. The polarized 
synchrotron radiation of the Galactic non-thermal radiation also serves to probe 
the ISM, including space between the emitting regions and the solar system. Radio 
propagation measurements provide unique information on turbulence in the ISM 
as well as the mean plasma properties such as density and magnetic field strength. 
Radio propagation observations can provide input to the major contemporary 
questions on the nature of ISM turbulence, such as its dissipation mechanisms 
and the processes responsible for generating the turbulence on large spatial scales. 
Measurements of the large scale Galactic magnetic field via Faraday rotation pro- 
vide unique observational input to theories of the generation of the Galactic field. 



Keywords interstellar matter-Milky Way, 98.38.-j ■ plasmas-astrophysical, 
95.30.Qd ■ turbulence-space plasma, 94.05.Lk 



1 Introduction 

The purpose of this article is to discuss how radio propagation measurements 
provide diagnostics of the interstellar medium (ISM) . By radio propagation mea- 
surements, we mean those in which a radio astronomical observable (such as the 
interferometric visibility, or the polarization position angle) has been modified by 
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a medium between the source of the radio waves and the radio telescope. In this 
paper, we will be interested in plasma media. These measurements provide rather 
direct information on the ionized gas density (strictly speaking, the electron den- 
sity) , the interstellar magnetic field, and (sometimes) indirect information on flow 
velocities in the interstellar medium. 

In addition to information on the mean plasma properties of the interstellar 
medium such as (n e ) and (B), these propagation observations yield information 
on turbulent fluctuations in the interstellar plasma. In fact, it can be argued that 
the information on interstellar plasma turbulence is the most unique contribution 
of this type of observation to studies of the ISM. 

This paper is intended, in part, to serve a tutorial and review function. How- 
ever, there have been numerous reviews in the past on the probing of the inter- 
stellar medium by radio propagation measurements, and the impli cations of those 
measurements for th e astrophysics of the ISM (see, in particular. lUscinski 1119771 
iRickett I Il977l Il990l) . There is no point in repeating the material already pub- 
lished in those papers. In the present article, we will make detailed reference to 
those papers to make a number of important points. At the same time, we will 
stress remaining, open questions about the interstellar plasma and its turbulence. 
In some cases, those questions have been actively discussed for many years. We 
will also clearly point out and discuss those topics in which radio propagation 
measurements provide crucial data for some of the issues of greatest importance 
in contemporary astrophysics. 

An underlying theme of this paper will be the conceptual unity of plasma pro- 
cesses that occur in the interstellar medium, the solar corona, the interplanetary 
medium, and finally, experiments in plasma physics laboratories. In the last decade 
or two, plasma physics laboratory experiments have succeeded in illuminating 
processes which also occur in astrophysical plasmas. These experiments deal with 
processes which are at the basis of plasma astrophysics. A partial list of the experi- 
ments which are contributing a new dimension to plasma astrophysics are measure- 
ment of Faraday rotation in laboratory plas mas, and its use in d i agnosing the basi c 
properties and processes in those plasmas ( Brower et al l 2002t Ding et al 1120031) . 



observation of the nonlinear interaction of Alfven waves ( Carter et al 1120061) , and a 
number of experimental efforts to investigate the nature of magnetic field reconnec- 
tion, a core process in astr ophysics (e.g. iBrown et al 1120021 ; IZweibel and Yamada I 
l2009t lYamada et al 1120101) . The unity of plasma physics and plasma astrophysics 
is exemplified by the interesting fact that the same radio propagation techniques, 
with the same radio telescopes, are, or can be, used to study the plasma physics 
of the interstellar medium, the corona, and the solar wind. 



1.1 The Fundamentals: (1) Phases of the ISM 

As has been noted for decades, the interstellar medium exists in a number of 
"phases" of different temperature, density, and ionization state. These different 
properties mean that fundamental plasma parameters such as the ion gyroradius, 
ion cyclotron frequency, plasma /?, and Debye length differ from one phase to 
another. The different phases and their plasma parameters were discussed in, and 
summarized Table 1 of ISpangler I (|200lh . For ease of reference, this information is 
reproduced here as Table 1. An interesting recent contribution to the discussion 
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Table 1 Main Phases of the Interstellar Medium 


Astronomical Name Density (cm -3 ) 


Temperature (K) 



Molecular Cloud 


> 10 4 


< 70 


Cold Neutral Medium (CNM) 


10 - 100 


~ 100 


HII regions 


5 - 10 


8000 


Warm Neutral and Ionized Mediua (WNM & WIM) 


0.1 - 0.5 


8000 


Coronal (HIM) 


io- 3 


10 6 



of phases of the I SM has been the advocacy of Heiles for a "fifth phase" , which is 
partially ionized ( Heiles 1120111) . We concur with Heiles that even a relatively small 



degree of ionization can profoundly alter the dynamics of a largely neutral gas. 
Many laboratory plasmas have a lower degree of ionization than some of the neutral 
phases of the ISM listed in Table 1. However, it seems tha t this " fi fth ph ase" is 
largely the same as the WNM. The new development in iHeiles I ( 20111 ) is the 



recognition of the role played by the partial ionization. As pointed out there, all of 
the phases are probably plasmas, since even the cold, predominantly neut ral phases 



have sufficient ionization to satisfy the traditional criteria of plasmas ( Spangler I 
1200 if) . 

As has long been noted, the pressures of the less dense phases of the ISM are, 
very roughly, compa rable at a value of 1.0 x 10~ 13 — l.Ox 10~ 12 dynes/cm 2 ( Ferriere I 



11998c iTielens 1120051) . By this standard, molecular clouds and HII regions are over- 
pressured. In the case of molecular clouds, the gravitational potential contributes 
to confinement of the gas. HII regions are overpressured, expanding entities. The 
other phases, Cold Neutral Medium (CNM), Warm Neutral Medium (WNM), 
Warm Ionized Medium (WIM) , and Coronal Phase or Hot Ionized Medium (HIM) 
have roughly comparable pressures and may, in fact, be in pressure equilibrium. 
In addition, the magnetic pressure of the interstellar medium is comparable to 
the aforementioned gas pressures, with a value of ~ 1.4 x 10~ 12 dynes/cm for 
Bism = 6 x 10 -6 G. Finally, the pressure corresponding to the energy density of 
the Galactic cosmic rays is also similar, ~ 1.0 x 10~ 12 dynes/cm 2 , suggesting equili- 
bration between the forms in whic h the ISM can s tore energy. This whole situation 
is summarized in the textbook bv lTielens I (|2005l) . where a value of ~ 0.5 x 10~ 12 
dynes/cm 2 is quoted for the gas phases of the ISM (CNM, WNM, WIM, and HIM), 
and a pressure of ~ 1.0 x 10~ 12 dynes/cm 2 is assigned to both the magnetic and 
cosmic ray pressure. 

The pressures and other properties of the various phases of the ISM, as well 
as the press ures of th e inter st ellar magnetic fi eld and cosmic rays were considered 
in detail bv IFerriere I ( 1998(1 . IFerriere I (J19981 ) also estimated how these pressures 



would change with Galactocentric radius and altitude above the Galactic plane . 
For the Galactic location of the Sun, and in the Galactic plane, IFerriere I ( 19981) 
estimates (see Figure 3 of that paper) a pressure of 6.0 x 10~ 13 dynes/cm 2 for the 
gas phase, and ~ l.Ox 10~ 12 dynes/cm 2 for both the magnetic and cosmic ray 
pressures. 

The above considerations are relevant to the scope of this paper. Our view of 
the interstellar medium considers it as a dynamic plasma. Magnetohydrodynamics 
(MHD) describes the dynamics as an interaction of the gas and the magnetic field 
via pressure terms as discussed above and magnetic tension forces. In addition, 
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plasmas interact with energetic particles such as the cosmic rays through resonant 
interactions with plasma turbulence. 

Among the important, recent developments in this field has been continued 
progress in specifying the s trength of the interste llar magnetic field T$ism i an d its 
dependence on gas density ( Crutcher et al II2010T ) . A plot of magnetic field strength 



(largely deduced through Zeeman effect measurements) versus gas density shows 
considerable scatter, and a trend towards la rger values only for densities greater 
than about 300 cm -3 . ICrutcher et al I (|2010l) also infer a median magnitude of the 
interstellar magnetic field of 6 ± 1 /iG. This is in agreement with equipartition 
estimates of the total magnetic field strength from synchrotron emission, and a 
factor of about three higher than the regular magnetic field component. At the 
present, there is no observational evidence for a change in the magnitude of T$ism 
between different parts of the low density ISM. 



1.2 The Fundamentals: (2) Radio Wave Propagation through the ISM 

This paper will concentrate on two radio propagation effects, acting on small 
(10 2 — 10 3 km) scales and large (pc) scales, respectively. These are angular broad- 
ening of a compact or pointlike source due to density turbulence in the interstellar 
medium, and Faraday rotation of linearly polarized radio waves from a radio source 
embedded in the ISM, or from outside the Galaxy. We also briefly allude to other 
radio scintillation phenomena caused by small scale turbulence. In the latter topic, 
we include the signature of ISM turbulence in gradients of the polarization vector 
of synchrotron radiation. 



1.2.1 Angular broadening of compact sources 

The basic physics content of radio wave propagation through the ISM is to be 
found in the expression for the refractive index of radio waves in a plasma. This 
is discussed in the proceedings of previous meetings of the International Space 
Science Institute, i.e. Spangler (2001) and Spangler (2009). As discussed there, 
the refractive index depends on the plasma density and magnetic field. For radio 
wave propagation in the ISM, the magnetic field dependence is determined by 
the component of the magnetic field in the direction of wave propagation. The 
modification of the radio refractive index by the magnetic field is much smaller 
than the modification due to the plasma density. This is responsible for the well- 
known feature that radio propagation measurements primarily diagnose the plasma 
density of the ISM, with only a higher order contribution due to the Galactic 
magnetic field. 

Turbulent fluctuations in the plasma density and magnetic field cause stochas- 
tic spatial and temporal fluctuations in the refractive index in the ISM. As a 
result, propagation through such a medium induces all manner of fluctuations in 
the received radio wave field ( Uscinski Ill977|) . The theory of how fluctuations in the 



refractive index generate corresponding fluctuations in properties of the wave elec- 
tri c field (va r ious n -point correlations of the electric field) is generally attributed 
to iTatarski I (|196ll) . An excellent illustration of the effects of wave propagation 
through a random medium is given by dynamic spectra of pulsars, an example of 
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Fig. 1 A dynamic spectrum of the pulsar PSR1133+16. Wave propagation through the 
stochastic interstellar medium produces constructive and destructive interference at different 
spatial locations (here portrayed as time) and frequencies. Observations and figure provided 
by James Cordes, Cornell University. 



which is shown in Figure[][fj The spectrum of a pulsar is a measured as a function of 
time, and the set of spectra combined as shown in Figure [1] In the absence of the 
turbulent interstellar medium, the flux density of the pulsar would be constant 
over the frequency range shown. The gray scale indicates the brightness of the 
pulsar, with dark shaded regions being bright. The variation in brightness is due 
to scattered radio waves alternatively constructively and destructively interfering 
at different frequencies and times. A discussion of pu lsar dynamic spectra and 

(1986|). The specific observations 



the information they contain is given in Cordes 



shown in Figure 1 are discussed in lLazio et al I ( 



A major goal of the theory of wave propagation in a random medium is to 
relate, via an integral transform, the radio astronomical measurement as a function 
of the independent variable, to the density power spectral density as a function 
of wavenumber. Examples of m ajor contribu t ions to this literature are lUscinski I 
(|l977T ). lLee and Jokipiil l|l975l) . and lRickettl (|l977lfl99oh . An illustration of the 
va rious type s of ob servable stochastic propagation phenomena is given in Figure 1 
of lSpangler I (120091) . 



We thank James Cordes of Cornell University for providing this graph. 
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1.2.2 Depolarization and Faraday Rotation of Synchrotron radiation 

Interstellar radio propagation measurements also provide information on the basic 
plasma state of the ISM, such as the plasma density, the vector magnetic field, as 
well as how these fields vary with position in the Galaxy. 

Since variations in the magnetic field vector along the line of sight and/or 
within the angular size of a telescope beam will partially depolarize linearly polar- 
ized synchrotron emission, the observed degree of polarization traces the ra tio of 
large -scale (regular) magnetic field strength to total magnetic field strength ( Beckl 



120011) . However, due to small-scale variations in this ratio due to local structure 
(supernova remnants, variable turbulence parameters), this method is mostly uti- 
lized on kpc-scales in external galaxies. 

Parsec-size scales in the magnetized ISM are typically probed using Faraday 
rotation. The Faraday rotation measure RM is directly proportional to the path 
integral along the line of sight (los) of the electron density n e and line-of-sight 
component of the magnetic field B« : 

Measurements of RM therefore provide nearly unique information on the magnetic 
field in the tenuous, ionized component of the ISM. 

Traditionally, Faraday rotation is measured from the rotation of linear po- 
larization angle 8 as a function of wavelength 6 = 9q + RM\ 2 , where 9o is the 
intrinsic polarization angle at emission of the synchrotron radiation. An illustra- 
tion of these ideas is given in Figure [2 which shows the Faraday rotation measure 
RM along several lines of sight to extragalactic radio sources in the Galactic plane 
in Cygnus. The Faraday rotation of the synchrotron radiation emitted by these 
sources is dominated by the Milky Way. The large differences in magnitude, and 
even sign, of .RM between closely-spaced lines of sight in Figure[2]are an indicator 
of the role of young, luminous stars in this region, as they produce structure in 
the ISM via stellar winds and supernovae, and ionize the gas. 

However, if synchrotron emission and the Faraday-rotating medium are mixed 
or alternating along the line of sight, the simple linearity of polarization angle 
change with A 2 is no longer valid. This may be the case in the majority of Faraday 
rotation measurements. Faraday rotation measurements of the diffuse synchrotron 
emission in galaxies is the most obvious example, but extragalactic sources may 
also have several intrinsic RM components. In this case, every RM component 
i along the line of sight - now called Faraday depth 4>i to indicate that it only 
probes Faraday rotation along a part of the line of sight - adds its own polarization 
angle rotation of <f>i\ , resulting in a non- linear polarization angle change with A . 
However, this opens the possibility of a Fourier transform, with A 2 as one of the 
conjugate variables, in order to disentangle the various <f>i components in the total 
observed signal. This metho d is called Rotation Measure synthesis ( Bur nl [1966; 
iBrentiens fc de Bruvnll2005h . 

This Fourier transform relation between observed polarization vector as a func- 
tion of wavelength squared P(A 2 ) and the Faraday dispersion function (or Faraday 
spectrum) is expressed as 



Plasma Diagnostics of the Interstellar Medium with Radio Astronomy 



Intensity [R] 



o 

□ 




74 72 

galoctic long 

Fig. 2 Highly variable Faraday Rotation Measures in the Galactic plane in Cygnus, in the 
vicinity of the Cygnus OBI star formation region. The sizes of the plotted circles are pro- 
portional to the magnitude of RM, specifically proportional to the logarithm of the absolute 
magnitude of RM. Filled circles correspond to positive RM, open circles are negative RM. The 
largest value of (absolute magnitude) of RM is 732 rad/m 2 , and the smallest is 7.5 rad/m 2 . 
The gray scale represents the intensity of Ha emission in this region; the Ha data may show 
the upper half of the plasma shell associat ed with the Cygnus OBI associ ation. Further discus- 
sion of these data is given in Whiting ct al (2009). Figure taken from Whiting ct al (2009), 
reproduced by permission of the American Astronomical Society. 
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However, since integration over wavelengths from — oo to +00 is not possible 
by definition, in practice these equations include a window function W{\ 2 ) which 
is non-zero where there is wavelength coverage and zero elsewhere: 
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This introduces the Rotation Measure Spread Function (RMSF) R(4>), which 
describes sidelobes in the Faraday depth signal due to imperfect wavelength cov- 
erage, very analogous to the dirty beam (or point spread function) in radio inter- 
ferometry due to imperfect coverage of the (u, v) plane. 
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Fig. 3 Polarized intensity (top) and polarization angle (2nd from top) as a function of wave- 
length squared for synchrotron emission from an unnamed cxtragalactic point source. The 
Faraday spectrum (3rd from top) as a function of Faraday depth shows two peaks, indicat- 
ing two separate Faraday depth components along this line of sight. Credit: Schnitzeler et al, 
A & A 494, 611, 2009, reproduced with permission ©ESO. 



A nice illust r ation of this effect is given in Figure [3l which is taken from 
ISchnitzeler et all ( 20091) . The Figure shows polarized synchrotron intensity and 
polarization angle as a function of wavelength squared of an extragalactic point 
source. Variations in polarized intensity and non-linearity of polarization angle 
with A 2 suggest that synchrotron emission and Faraday rotation are partially 
mixed. Indeed, the Faraday spectrum in the Figure shows two Faraday depth peaks 
belonging to two different synchrotron emitting regions with different amounts of 
Faraday rotation along the line of sight. The two Faraday depth components might 
both be due to the Galactic interstellar medium, or one might belong to a magneto- 
ionized medium in or around the extragalactic point source. Without additional 
information it is not possible to distinguish between these two possibilities. 
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2 Different Astrophysical Media, Common Phenomena 

We point out the remarkable fact that a number of "astrophysical plasmas" , i.e. 
distinct and very different media, have roughly similar radio propagation effects. 

— The Faraday rotation measure (RM) due to the Earth's ionosphere is typically 
in the range 0.5 — 3 rad/m . This is of the order of the RM due to the solar 
corona for lines of sight that pass at a heliocentric distance ~ 10 Rq, and a 
factor of a few smaller than the st andard dev iatio n in R M of the interstellar 
medium at high Galactic latitudes iMao et all fe.g. I2010T) . 

— Very Long Baseline Interferometers operating at frequencies of 1 - 5 GHz mea- 
sure similar effects, and of similar magnitude, when observing extragalactic ra- 
dio sources through the inner sol ar wind at heliocentric distances of 10 - 30 Rp, 
or through the Galactic plane (Spangler a nd Sakurai I 1995t Spangler et al I 



l2002HSpangler and Cordes II 19881 iFev et allll99ll;ISpangler and Cordes II 19981) 



The fact that similar observational effects are measured for very different media 
means that scientists who study the ISM should be in dialog with heliospheric 
scientists who are studying related scientific questions. 

In the case of Faraday rotation, there is also a connection with laboratory 
plasma research. Faraday rotation is used as a diagnostic of fusion plasmas 
(|Brower et"aTll2002l; Iping et al 1120031) . This offers ISM astronomers the possi- 



bility of laboratory "ground truth" for some of our diagnostics. 



3 Well-Established Results from Radio Propagation Studies 

There are several results from radio propagation studies that are so well estab- 
lished, and confirmed by different investigators, that they now have the status of 
basic properties of the ISM that should be explained by theories. 



3.1 "The Big Power Law in the Sky" 

One of the best-known results from radio propagation studies is that the spatial 
power spectrum of density fluctuations follows a Kolmogorov spectrum over at 
least 5, and perha ps 10 decades (lArmstrong et al Ill995l) . It should be emphasized 
that the result of lArmstrong et al I (J1995I) pertains to the Warm Ionized Medium 
(WIM) component of the ISM. Analogous results hav e been reported for diag- 
nostics of the neutral gas (e.g. IChepurnov et al 1 120101) . The neutral gas resides 
in the WNM phase, which is spatially distinct from the WIM. As noted above, 
(Section 1.1) the degree of ionization in the WNM, while low, is not zero, and it 
m ay be considered a s a we akly-ionized plasma. Whether the fluctuations studied 



bv IChepurnov et al I (2010) are driven by plasma dynamics due to the minority 



ion species remains unknown. 

It is intriguing that the solar wind density power spectrum is the same as 
that observed for the ISM (at least for the slow solar wind), although t he inertia! 
subrange is smaller, being perhaps 3 orders of magnitude (Bruno and Carbone 
120051) . In the case of the ISM, it is perhaps not sufficiently appreciated that there 
must be a corresponding power law spectrum of magnetic field fluctuations to 
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explain the highly diffusive transport properties of cosmic rays w ith a very wide 
range of energies. This point has been made and emphasized bv Ijokipii I ( 19771 
Il988h . 



3.2 Spatial Variation in the Intensity of Interstellar Turbulence 

A measure of the intensity of turbulence is the parameter Cfj, which is the nor- 
malization constant of the density power spectrum. That is, if the spatial power 
spectrum of density fluctuations is P n (q) where q is the spatial wavenumber, 

Pn(q) = C 2 N q- a (7) 

and a is the spectral index of the power spectrum. This form of the power spec- 
trum is the simplest, in which the power spectral density depends only on the 
magnitude of q, and is thus isotropic. It is known that magnetohydrodynamic 
(MHD) turbulence is in fact anisotropic, with the large scale magnetic field deter- 
mining the preferred direction (See comments in Section 3.4 below). In this case, 
the power spectral density P n (q) depends separately on q± and qu, the compo- 
nents of q pe rpendicular an d parallel, respectively, to the large scale field (this is 
discussed in ISpangler Ill999l with a guide to the relevant literature). For most of 
this paper, we adopt EquationQas a convenient approximation. 

The parameter C^r is directly related to the variance of the density fluctua- 
tions. It has been long realized that C% varie s drastically from one part of the 
inter stellar medium to another ( Rickett Ill977t ICordes. Weisberg. and Boriakoff I 



ll985T) . In some cases, it is clear that lines of sight with large C% traverse HII 
regions, or other regions with higher than normal plasma density. However, there 
appear to other lines of sight where no such obvious region of enhanced density 
exists. It remains uncle ar whether some of this variation could be due to true 
turbulent intermittency ( Spangler and Cordes Ill998t ISpangler Ill999j) . 



3.3 The Galactic-Scale Magnetic Field 

Large-scale Faraday rotation surveys of the sky show substantial organization of 
the RM in different parts of the sky, in the sense that the magnitude and sign of 
RM are correlated over significant parts of the sky. This is interpreted as evidence 
of a Galactic-scale magnetic field. Most likely, the large-scale magnetic field in the 
Milky Way generally follows the sp iral arms, as is ubiquitously seen in synchrotron 
observations of external galaxies ( Beckl 120011) . However, there is some evidence 



for local deviations from the spiral structu re (IBrown et al. l2007t iRae fc Brownl 



l2010t IVan Eck et"alll201ll) . similar to M51 (jPatrikeev et al.l 120061) . One reversal 
in the large-scale magn etic field direction just in side the Solar circle has been 
known for decades (e.g. iThomson fc Nelson! Il980l) . although it remains unclear 
why these large-scale reversals are not observed in external spirals. There is still 
much controversy about the number and location of any other larg e -scale field 



reversals in the Galactic disk. For an extensive review, see iHaverkornl ( 20131) . 
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3.4 Anisotropy of Turbulence 

A major result from the theory of magnetohydrodynamic (MHD) turbulence, con- 
firmed by observations of turbulence in the solar corona and solar wind, is that 
turbulence is anisotropic, in the sense of turbulent irregularities being stretched 
out al ong the large scale magnetic field. This result was obtained by IStrauss I 
( 1976J) for the case of irregularities in fusion plasmas, but the arguments presented 



by Strauss are also valid in the case o f astrophysical plasmas. This prop erty is also 
present in interstellar turbulence fsee lBrisken et a l 2010; Rickett 11201 lL for recent 
discussions of the more pronounced cases). 



3.5 The Dissipation Range of Turbulence 

Important progress has been made in the past decade in our understanding of 
the dissipation of plasma turbulence. This progress has been possible through 
improved m easurements of solar wind turbulence, as wel l as novel theoretical de- 
velopments dHowes et al Il2008t TAlexandrova et al I [20091: iHowes et~aT1l201lll201lt 
ISahraoui et al II2012T) . These investigations have identified spatial scales on which 



dissipation occurs, and advanced suggestions for the mechanisms by which this 
occurs. It is now clear that a break in the power spectrum of magnetic field fluc- 
tuations in the solar wind occurs on scales comparable to, and smaller than the 
ion inertial length li, 

>>•% < 8 » 

where Va is the Alfven speed and J? ? - is t he ion (proton) cyclotron frequency 
( Howes et al I l2011t lAlexandrova et al 1 120091) . There remains active discussion in 



the community as to whether the dissipation is due to Landau damping of highly 
oblique Alfven waves (the assumption being that turbulent fluctuat ions on these 
scale s have damp ing properties similar to linear plasma wave modes ( Howes et al I 
120081 l201l|J201l|)), dissipation of other, higher frequen cy modes propagatin g at 
large angles with respect to the mean magnetic fields (Sahraoui et al 2012), or 
damping by other modes on electron inertial scales (JAlexandrova et al 1 120091) . In 
keeping with the philosophy of this paper, we assume that these results are of 
great importance of our understanding of the ISM as well. 

As will be discussed in Section 4, there is now good evidence for the beginning of 
the dissipation range in interstellar turbulence at the ion inertial length, although 
there is also evidence for differences between solar wind and interstellar turbulence 
(see Section 4.5 below). 



3.6 The Outer Scale of Interstellar Turbulence 

The outer scale to magnetized turbulence in the Warm Ionized Medium (WIM) 
phase of the ISM is measured from structure functions of RM to be a few parsecs 



in the spira l arms, but up to ~ 100 pc in the interarm regions ( Haverkorn et al I 
2006, 2008) . This is in agreement with estimates of the outer scale of turbulence 
averaged over large parts in the sky ( mostly towa rds the Galactic halo) of or - 
der 100 pc by iLazarvan fc Shutenkovl (e.g. Il990h : lOhno fc Shibatal (e.g. I1993J) : 
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IChepurnov and Lazarian I fe.g. l201dh . These scales are similar to the fin al sizes of 
supernova remnants, suggesting (combined with energy arguments, see iMacLowl 
(2004)) that these are the dominant sources of turbulence in the WIM. The smaller 
outer scale found in spiral arms may be due to the abundance of H n regions 
( Minter fc Spangler Hl996tlHaverkorn et al II20041 ). Small out er scales of a few par - 



secs have also been found in the highly polarized Fan region j lacobelli et al.ll2013ri 
or from anisotropies in TeV cosmic ray distributions ( Malkov et al.ll2010l) . 



4 Turbulent Microscales in the Interstellar Medium 

4.1 Definition of Turbulent Microscales 

In this section, we discuss the ways we can measure turbulence on very small scales 
in the interstellar medium. By very small scales, we mean those in the dissipation 
range. Knowledge of this part of the turbulent cascade is very important because it 
contains information on the way in which energy is taken from the large scales, on 
order of parsecs or more, and transferred to other forms, presumably heat energy 
of the interstellar gas. 

Our interest in this section will be particularly focused on two phases of the 
ISM, the WIM and HII regions surrounding young stars. The WIM is of interest 
because it appears to be the best-diagnosed phase of the ISM, in the sense of a 
plasma physics measurement, primarily as a result of investigations done at the 
University of Wisconsin using Fabry Pe rot interfe rometers as very high spectral 
resolution spectrometers (see Irlaffner et alll2003l for a recent summary of this 
work.). 



4.2 Turbulent Microscales in the Solar Wind 

Once again, we use the solar wind, with its extensive and often sophisticated, 
in-situ measurements and substantial body of theoretical results as a model for 
interstellar turbulence. Spacecraft instruments provide in-situ measurements of 
virtually all plasma parameters of interest, and provide the best data set for dis- 
cussions of MHD turbulence. Measurements of solar wind turbulence at a helio- 
centric distance of 1 a.u. (the bulk of spacecraft measurements) show a power-law 
power spectrum of magnetic field fluctuations with a single spectral index that ex- 
tends from an outer scale with a size of one to a few solar radii, do wn to an inner 
scale of a few thousand kilometers (e.g. iBruno and Carbone 1120051) . A number of 
investigations have shown that this scale corresponds to the ion inertial length 1$ 
defined in Equation [8] 



R adio propagation observations through the corona and in ner solar wind (IColes and Harmon I 

Il989 : Harmo n and Coles 112005c ISpangler and Sakurai Ill995l) also show strong ev- 
idence for an enhancement in the power spectral density of plasma density fluc- 
tuations on the scale of the ion inertial length. This bulge on the approximate 
scale of the ion inertial length can also be seen in power spectra from in-si tu mea- 



surements of plasma density in the solar wind at 1 a.u. (jChen et al 2012). These 
observations, both from direct, in-situ measurements as well as radio propagation 
observations, are interpreted as evidence that the fluctuations in the dissipation 
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range have properties of obliquely-propagating, kinetic Alfven waves, since such 
waves become more compressive at the ion inertial scale (Harmon 1989). In fact, 
the kinetic Alfven nature of the fluctuations on the dissipation scale is the basis of 
the model of turbul ence in the dissipation range advanced by Howes and colleagues 
( Howes et alll201ll ). So, the situation for the solar corona and solar wind seems 



quite consistent as regards both measurements of plasma density fluctuations and 
theoretical understanding of the entire turbulent cascade. The prominence of this 
directly-detected bulge appears to be less pronounced than that retrieved from 
radio propagation measurements in the corona and inner solar wind. This may in- 
dicate that the kinetic Alfven wave component decays with increasing heliocentric 
distance in the solar wind. This would hardly be surprising , since many proper- 
ties o f the solar wind change with heliocentric distance ( e.g. iBruno and Carbone I 
l2005h . 



4.3 How We Measure Microscales in Interstellar Turbulence 

The ion inertial length in the WIM ph ase of the interstellar mediu m is of order 
hundred to a few hundred kilometers ( Soangler and Gwinn I ll99(J) . At first, it 



seems amazing that any kind of astronomical measurement, made on a medium 
with an extent of kiloparsecs, could diagnose fluctuations on such a scale. Radio 
propagation measurements make this possible. 

The subsequent discussion in this section will concentrate on one of the scin- 
tillation phenomena mentioned in Section 1.2, angular broadening. A point source 
of radio waves viewed through a turbulent medium will appear as a blurred, fuzzy 
object. Essentially the same phenomenon is encountered at optical wavelengths 
in the form of seeing disks of stars. A radio source that is blurred by interstellar 
turbulence will have a measured brightness distribution I(x, y) which is more ex- 
tended than the intrinsic image of the source. This brightness distribution contains 
information on the intensity and spatial power spectrum of the density fluctua- 
tions. Here I is the intensity of the radiation, which is a function of two angular 
coordinates on the sky, x and y, customarily Right Ascension and Declination. 
The brightness distribution is related to the observable quantity which is directly 
measured by the interferometer. The way a radio interferometer makes an image of 
a radio source is to measure the complex visibility function V(u, v), which is directly 
related to the correl ation between the radi o wave electric field at two antennas 
of an interferometer (JThompson et al Ill986l ). The complex visibility function has 
units of Janskys (radiative flux), and is a function of the arguments u and v, the 
east-west and north-south components of the interferometer baseline, normalized 
by the wavelength of observation. A two dimen sional Fourier transform relates the 
complex visibility function V(u,v) and I(x,y) (JThompson et allll986h . 

Observations have shown that the turbulent irregularities in the WIM and 
HII regions around OB associations, like those in the corona and solar wind, are 
anisotropic in the sense that is theoretically e xpected. A summary of the observa- 
tional evidence as of 1999 is given in Spaneler (1999). For the present purposes, 



we will employ the approximation of isotropic irregularities . The arguments given 
here can be generalized to the case of anisotropic scattering ( Spangler and Cordes I 



14 Marijkc Haverkorn, Steven R. Spangler 



ll988T) . In the case of isotropic scattering, 



V(u,v) = V(\/u 2 +v 2 ) =V(r) (9) 

where r is the (dimensional) interferometer baseline length, projected on the plane 
of the sky. In this case, the complex visibility function of a point source viewed 
through a turbulent medium is a real function, ( Cordes. Weisberg. and Boriakoff I 
ll985l:lSpangler and Gwhm"lll990l) 

V(r) = S e~* D * W (10) 

where So is the total flux density of the source, and D^,(r) is the phase structure func- 
tion, which contains information on the intensity and spatial power spectrum of the 
density fluctuations. For the case of scattering by a homogeneous, turbulent slab 



of th i ckness L, the phase structu re function is ( Cordes. Weisberg. and Boriakoff I 
Il985t ISpangler and Gwhm"lll990h 



D^r) = 8n'rtX"L / dqq [1 - J (qr)] P n (q) (11) 

Jo 

The functions and variables in Equation [11] are as follows. The classical electron 
radius is given by r e , A is the wavelength of observation, L is the thickness, or extent 
along the line of sight of the turbulent plasma, q is the magnitude of the turbulent 
wavenumber, Jq(x) is a Bessel function of the first kind of order 0, and P n (q) is the 
spatial power spectrum of the density fluctuations. One expects the power spectra 
of all plasma parameters to be modified at wavenumbers corresponding to the 
reciprocal of the ion inertial length, ion gyroradius, or similar plasma mic r oscale 
on which dissipation begins to become important. ISpangler and Gwinn I ( 19901) 



adopted the following simple model in which EquationQis modified by having the 
power spectrum truncated on wavenumbers larger than a dissipation wavenumber 
Qo, 

Pn(q) = C 2 N q- a e-«l'» (12) 

Although Equation[12]is highly simplified, and was adopted bv lSpangler and Gwinn I 
(1990) for analytic convenience, the power spectrum of magnetic field fluctuations 



in the solar wind at 1 AU is truncated by an exponential function (JAlexandrova et al I 
20091 see also presentation by A lexandrova et al in this m eeting). An important 
difference between th e result of Alexandr ova et al I ( 20091) and the analysis pre- 



sented below is that lAlexandrova et al I ( 2009J) found exponential truncation of 



the power spectrum on electron rather than ion scales (i.e. electron inertial length, 
or electron gyroradius). 

Substitution of Equation [12] into Equation [TT] and change of variables from 
q — > y = qr gives the following expression 

/>oo 

D^r) = 87v 2 r 2 e \ 2 (C 2 N L)r a - 2 dy [1 - My^y-^ 1 ^-^ (13) 

Jo 

where Q = j-, where l& is the dissipation scale, J<j ~ ~. The quantity C^L is 
termed the scattering measure, and roughly determines the magnitude of angular 
broadening. Turbulent plasmas that have large Cjy, large L, or both, will produce 
heavy angular broadening. 
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Fig. 4 Plot of the function I(y,Q) for a Kolmogorov spectrum (a = 11/3) and Q = oo. The 
integrand tells what wavenumbers in the turbulent spectrum dominate the measurement on a 
given interferometer baseline. This occurs for y = qr > 1, which corresponds to irregularities 
with sizes of order the interferometer baseline. 



Equation [T3] yields important insight on remote sensing diagnosis of inter- 
stellar (and heliospheric) turbulence. Let us start with the case Q — > oo, which 
corresponds to an infinitely small dissipation scale. In this case, the spectrum is 
power law for all wavenumbers larger than that corresponding to the outer scale. 
In this case, the integral in Equation 1131 is a number which depends only on the 
index a. For the Kolmogorov spectrum (a = 11/3) the value of the integral is 
1.117, and the structure function D<p(r) oc r 5 ' 3 . 

Equation [13] also illustrates one of the most intriguing aspects of radio wave 
propagation, and demonstrates why radio astronomical measurements can con- 
tribute much to a discussion of plasma microscales in the ISM. Since Equation [13] 
is an integral over wavenumber (this is explicit in Equation lll[) , the integrand 
shows which wavenumbers dominate the measurement. The integrand in Equa- 
tion [T3] which we note by the function I(y,Q), is shown in Figure [3] for the case 
I(y, oo) and a = 11/3 (Kolmogorov spectrum). 

The function I(y, oo) is monotonically decreasing with increasing y. At first, 
this would seem to indicate that the lowest wavenumbers in the spectrum dom- 
inate the measurement. However, that is not true for the case of a Kolmogorov 
spectrum. The measurement (D^) is determined by an integral over y. This is be- 
cause for y < 2, each progressively higher decade in y makes a larger contribution 
to the integral. The integral is dominated by contributions with y ~ 1 — 10 where 
there is an inflection in the function 7(y,oo). When an interferometer measures 
a broadened radio source, it is responding to irregularities with a wide range of 
wavenumbers. However, the dominant contribution to the visibility measurement 
is from irregularities with sizes comparable to the baseline length. This baseline 
length ranges from tens of kilometers in the case of the Very Large Array, to a few 
thousand kilometers in the case of Very Long Baseline Interferometers. Obviously, 
for this statement to be relevant, a measurement with a given interferometer must 
be affected or even dominated by p ropagation effect s. This point was made in the 
context of interstellar scattering in ISpangler I ( 1988|) . 
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Fig. 5 Theoretical structure functions of turbulence in the interstellar medium. The plotted 
lines correspond to theoretical structure functions D^r) given by Equation (10). Both adopt a 
Kolmogorov spectrum of density irregularities. The solid curve is for a Kolmogorov spectrum 
with no inner scale (l^ —> 0). The dashed line corresponds to a spectrum with an inner 
scale Id = 300 km. Sensitive, carefully calibrated interferometer measurements can distinguish 
between these two cases. 



Equation 1131 also shows how the presence of turbulent dissipation is manifest 
in radio propagation measurements. When Q is finite, corresponding to a finite 
dissipation scale, the e~ v ' < ^ term will depress the value of the integral. The value 
of D,p(r) at short baselines, where dissipation is pronounced, is less than a value 
extrapolated from larger values of r according to an r ' relation. In the dissipation 
range, Dj,(r) has a steeper dependence than r 5 ' 3 . To illustrate these points, Figure 
5 shows two structure functions, one without an inner scale and the other with an 
outer scale of 300 km. 



4.4 Observational Results on Turbulent Dissipation Scales in the Interstellar 
Medium 



These issues were discussed in Spangl er and Gwinn I ( 19901 ). who showed that ob 



servers who interpret their angular broadening data in terms of a spectral in- 
dex a would report values which depend on the baselines used in the measure- 
ment. Angular broadening measurements on short baselines would yield a value of 
a ~ 4.0, whereas measurements on long baselines in the inertial subrange would 
yield a ~ 3.67. This may be seen by reference to Figure 5. The inferred spectral 
in dex is determined by the slop e of D^r) versus r on a log- log plot, such as Figure 



5. ISpangler and Gwinn I ( 1990J) assembled the data on angular broadening mea- 



surements that were available at that time, and showed that a dependence of the 
inferred value of a on the interferometer base lines used did seem to be pre sent in 
the data. The result was shown in Figure 1 of ISpangler and Gwinn I (|l990h . From 
these data they found that there could be a break in the interstellar density power 
spectrum with an inner scale of 50 - 200 km. More importantly, they pointed out 
that a scale in this range was actually expected, if the inner scale corresponds to 
the ion inertial length li defined in Equation [8] as is the case for scattering in the 
corona and solar wind (Section 4.2). 

A more direct determination of an inner scale, made from comparison of mea- 
sured D^r) values with the theoretical expression in Equations (10) and (12), 
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was made by iMolnar et al I ( 1995T) . iMolnar et al I ( 1995 ) made and analy z ed an 



gular broadening measurements similar to those of lSpangler and Cordes I ( 19881) 



but of the radio source Cygnus X-l. It is viewed t h rough the HII region associated 
with the Cygnus 0B2 association. IMolnar et al I ( 19951) essentially made a fit of 



Equation (10) (but including the effect of anisotropy of scattering) to their data, 
and found a satisfactory model to be a Kolmogorov underlying spectrum and an 
inner scale of 300 kilometers. 

An additiona l , and particularly compelling result has been the recent report 
by iRickett et al I ( 2009T) . They studied the form of the broadening profile of pulses 



of the pulsar PSRJ1644-4559 late in the pulse. At this time, radiation is being 
rece ived from highly sca ttered rays that are probing very small scale irregulari- 
ties. IRickett et al I (|2009j) make the important point that the amount of radiation 
received late in the pulse is only consistent with a Kolmogorov spectrum that 
breaks at an inner scale as expressed by Equation 1121 or something similar. A 
power spectrum that remained Kolmogorov to infinitely high spatial wavenumber 
woul d cause more p ulse p ower to be observed late in the pulse than is actually 
seen. IRickett et al I ( 20091) use their data to extract a value for the inner scale of 



70 - 100 km. 

These three independent investigations using radio propagation data have 
therefore concluded that there is a spectral break in the density power spec- 
trum of the interstellar medium, and that this inner scale is consistent with the 
ion inertial length. It shoul d be emphasize d that the li nes of sight analyzed by 
ISpangler and Gwinn I (|l990l) . IMolnar et all (|l995l) . and IRickett et all (|2009l) all 
traversed HII regions. There is, as yet, no observational data that can determine 
the inner scale to the turbulence in the WIM. We do not know if such an inner 
scale would be at the ion inertial length. 



4.5 A Break or a Bulge? 

One interesting, and at the present preliminary feature emergent from these inves- 
tigations regards the transition to the dissipation range in interstellar turbulence. 
In HII region plasmas, the dissipation range appears to consist of a smooth steep- 
ening, without the bulge in the density power spectrum on the ion inertial length, 
as exists in the corona and solar wind (Section 6.1). Given the admittedly limited 
present information, it appears that the interstellar spectrum of density fluctu- 
ations has no compressive bulge at the inner scale. If confirmed by subsequent 
investigations, it could point to an important distinction between turbulence in 
the i nterstellar med i um a nd that in the solar corona and solar wind. The results 
from IRickett et al I (2009J) seem particularly compelling, because a bulge in the 



interstellar density power spectrum on the ion inertial scale would produce more 
pulse power at late arrival times than is actually seen (this point is clearly illus- 
trated in Figure 7 of IRickett et al I (J2009h '). 



If this bulge is missing in the interstellar density spectrum, what does it sig- 
nify? Does it imply that kinetic Alfven waves are not present in the turbulent field, 
or that the small scale irregularities in the interstellar medium do not evolve in 
a manner similar to kinetic Alfven waves? In that case, what is the nature of the 
fluctuations over such a large inertial subrange in the ISM? As menti oned in Sec- 
tion 4.2, the solar wind results may provide guidance; the results of IChen et al I 
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(2012J) indicate that the prominence of kinetic Alfven waves decreases with in- 



creasing heliocentric distance. Interstellar turbulence is comparatively much older 
in terms of the number of eddy turnover times, so it is certainly plausible that the 
kinetic Alfven wave component of ISM turbulence has dissipated. 

These ruminations need more extensive and more convincing observational 
demonstration. Fortunately, the instruments and observational techniques are op- 
erational and available. The instruments currently available for angular broad- 
ening meas urements are greatly improved over those used in the measurements 
cited aboyeJSpangler and Cordes Ill988t ISpangler and Gwinn Ill990t iMolnar et al I 
Il995t ISpangler and Cordes Ill998l) . Those investigations used Very Long Baseline 



Interferometers with much smaller bandwidths and correlator capability than are 
now available with the Very Long Baseline Array (VLBA) of the National Radio 
Astronomy Observatory (NRAO). In addition, the LOFAR low frequency radio 
telescope in Europe is now operational and has the c apability of making novel 
angular broadening measurements. Finally, the work of iRickett et al I (|2009|) also 
demonstrates the advances that have been made in pulsar measurements of the 
ISM, utilizing new, state-of-the-art pulsar processors on large single dish telescopes 
such as the Parkes antenna or the Green Bank Telescope of NRAO. Future inves- 
tigations with these powerful new instruments could illuminate the interesting 
question as to whether turbulent dissipation processes are the same in the solar 
wind and the plasma components of the interstellar medium. 



5 Turbulent mesoscales in the interstellar medium 

5.1 Rotation Measure Synthesis results 

The interpretation of three-dimensional Faraday depth cubes (i.e. intensity maps 
in spatial coordinates where Faraday depth is the third dimension) is anything but 
straightforward . 

In addition to the artefacts introduced by a non-Gaussian rotation measure 
spread function, as explained in Section 1.2.2, a number of other effects contribute 
to the difficulty of translating rotation measure cubes into physical properties of 
the interstellar medium. 

Firstly, in analogy to aperture synthesis, a limited range in wavelength squared 
causes a limited sensitivity to large-scale Faraday depth structures. In contrast 
with aperture synthesis, in RM synthesis this can lead to a situation where the 
maximum detectable scale is smaller than the Faraday depth resolution. Therefore, 
only Faraday-thin component^ and sharp gradients in Faraday depth, such as the 
edges of Faraday-thick components, will show up in a Faraday spectrum. The 
dependence on wavelength range of the Faraday depth resolution S<f>, maximum 
de tectable scale A4> max and ma ximum detectable Faraday depth 4>max are given 
bv lBrentiens &: de Bruvr] (|2005h as 



2 Faraday-thin component is denned as a gaseous medium observed at a wavelength where 
the ch ange of polarization angle through Faraday rotation is small. iBrcnticns &: de Bruvnl 
(2005) define Faraday-thin as </<A 2 << 1 and Faraday-thick as <j>\ 2 >> 1. A Faraday-thin 
component displays negligible internal Faraday depolarization. 
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Secondly, different Faraday depth features in a Faraday spectrum only contain 
information about the amount of their Faraday depth, but not necessarily about 
their distance. If along a line of sight Faraday depth increases monotonically, i.e. 
if no magnetic field reversals exist along the line of sight, then the distance order 
of Faraday depth components is known. However, in the general ISM, with many 
multi-scale magnetic field reversals, distance to Faraday components is usually 
unknown. Only in exceptional cases, if one has complementary rotation measures 
from a number of pulsars with known distances along similar lines of sight, or if 
the Faraday depth component has a counterpart with a known distance in another 
tracer, is it possible to estimate the distance to a Faraday component. 

Taking these caveats into account, a number of studies on RM synthesis of dif- 
fuse Ga lactic synchrotron emission have been done, which show mostly consistent 
results. iBrentiens I ( 20111) examines a ~ 4° x 7° field around the Perseus galaxy 



cluster, which mostly displays Galactic synchrotron emission, at a broad frequency 
range around 350 MHz. Synchrotron emitting components are detected at multiple 
Faraday depths between —50 and +100 rad m~ 2 , are Faraday thin and spatially 
thin (< 40 pc), and are well separated in Faraday depth space, suggesting that 
they are flanked by Faraday-rot ating-only parts of the ISM. 

The same effect is noticed by Iacobelli et al. (2013 ), who study the "rin g stuc- 
ture" in the Fan region (see e.g.Haverkorn e t alll2003t iBernardi et alH2009T ) in RM 



synthesis around 150 MHz. They also identify separate Faraday depth components, 
viz. the ring structure a nd a fo r egrou nd component which they associate with the 
Local Bubble. Similarly. IPizzol (2010j) notice in their Galactic foreground studies 



in the direction of the galaxy cluster Abell 2255 at multiple frequency bands from 
150 MHz to 1200 MHz three distinct ranges of Faraday depth with widely different 
morphologies. 

These early RM synthesis studies of Galactic diffuse synchrotron emission con- 
sistently conclude that the synchrotron emission is visible in detected in discrete, 
often Faraday thin, structures with widely different morphologies, interspersed 
with Faraday-rotating-only components. It is tempting to interpret these obser- 
vations as actual small-scale variability in synchrotron emission in the ISM, or 
discrete regions of excess emission. However, two other effects are at play as well. 
Synchrotron emission dominates in locations where B = B^, while Faraday rota- 
tion only depends on Bu. This may cause the observed apparent anti-correlation 
between synchrotron emission and Faraday rotation. Secondly, the insensitivity of 
the technique to large Faraday-thick (emitting and Faraday rotating) structures, 
which may mimic Faraday-thin emission components at the edges of the Faraday 
depth range, will play a role. 

Wavelet analysis can be successfully applied to Faraday depth cubes to recog- 
nize magnetic features such as turbulence or large-scale magnetic fiel d reversals 
in nearby spiral galaxies or the intracluster medium in galaxy clusters ( Beck et all 
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[2012J). Low frequency data (~ 100 MHz) are needed to provide the necessary Fara- 
day depth resolution, while broad frequency coverage (up to several GHz) is crucial 
to make broad Faraday structures detectable. In practice this requires combina- 
tion of broad-band data fr om various telescopes , such as the Global Magneto-Ionic 
Medium Survey f GMIMS. IWolleben et alll2009l) . 



5.2 Polarization gradients 



Linearly polarized intensity maps of diffuse synchrotron emission consistently show 
narrow one-dimensional structures of complete depolarization named depolarization 
canals (Haverkorn et al. 2000). Some of these depolarization canals are observa- 
tional artefacts due to missing short spacings in radio interferometric observa- 
tions, while other canals point to locations of sharp jumps in rotation measure, 
i .e. sudden changes in electron density and/or parall el magnetic field in the ISM 
( Shukurov fc Berkhuiisen 2003; H averkorn et al.l2004[ ). In addition, not all of these 
sudden changes in ISM conditions are visible as depolarization canals. 

The method of gradients in linear polarization was devised to ob tain a complete 
censu s of these locations of sudden change of conditions in the ISM ( Gaensler et al.l 
20111) . The vectorial polarization gradient is calculated from the Stokes parameters 
(Q, U) as 
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This gives a random-looking pattern of mostly one-dimensional locations of 
high polarization gradient, of which depolarization canals are a subset (see Fig. [6J. 
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Fig. 6 The map of gradient of linear polarization |VP| as denned in equation 1171 for a field 
around the Galactic plane. The inset panel shows a 0.9° X 0.9° close-up of the brightest structure 
in |VP|, where the direction of the gradient is plotted only for s trong (> 5 mjy beam 1,5 ) 
gradient amplitudes. Figure reproduced from lGaensler et al.l lt201lT) . 
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These polarization gradient filaments can be characterized by the moments 
of the polarization gradient distribution. Simulations of magnetohydrodynamic 
(MHD) turbulence show that the third and fourth order moments (skewness 7 
and kurtosis 0) increase monotonically with Mach number and depend on sonic 
Mach number. Comparison of the observed values 7 = 0.3 and /3 = 0.9 with 
simulated values for varying sonic Mach numbers indicates that the magnetic 
turbulence in the ISM (at least in the field given in Figure [6]) is mildly subsonic to 



transonic (jGaensler et al.ll201l() . This is in agreement with estimates of the sonic 



Mach number i n the warm ionized medium from emission measure distributions 
(JHill et alll2008h . 



Simulations also show that the filaments in polarization gradients can be 
caused by either inte racting shocks or random fluctuations in MHD turbulence 
( Burkhart et al II2012T) . These authors also introduce the genus method to charac- 



terize the polarization gradient maps. For subsonic turbulence as in the ISM, where 
magnetic field fluctuations dominate the polarization gradient topology, the topol- 
ogy is 'clumpy', as opposed to supersonic turbulence which shows a 'swiss cheese' 
topology. 



6 Mysteries of Interstellar Turbulence 

In the paper to this point, we have reviewed the remarkable amount of information 
revealed by radioastronomical studies on turbulence in the ISM. However, there 
remain a number of phenomena and effects that are not understood. In this section, 
we present what may be considered an agenda for future ISM radio propagation 
studies, that might clarify these issues. We discuss topics in which additional 
measurements may contribute in a major way to advances in our understanding 
of interstellar turbulence, or cases in which emerging observational results appear 
difficult to understand, given our current vision of the interstellar medium and 
the turbulence in it. These "mysteries" often involve input from the theory of 
plasma turbulence, and frequently rely on the latest results in this field. In other 
cases, they represent known observational results of long standing that have eluded 
explanation. 



6.1 The Existence of a Cascade in Interstellar Turbulence 

Is there really a cascade in interstellar turbulence from the outer scale of 4 parsecs 
(or 100 pc) to the dissipation scale of order 100 — 500 km? Even in the case of 
well-studied solar wind turbulence this issue is not entirely resolved. Observations 
of solar wind turbulence indicate that it is comprised of Alfven waves propagating 
in both directions with respect to the large scale interplanetary magnetic field, 
i.e. towards and away from the SurQ Within the context of the most commonly 
discussed theories of solar wind turbulence, these counterpropagating waves are 
necessary for the existence of nonlinearities that produce the turbulent cascade. 
In the case of the interstellar medium, we do not have observations sufficient 
to say if counterpropagating Alfven waves are present, and it is not clear how 



3 A more precise and technical statement would be that both the positive and negative 
"Elsasser Variables" are present in solar wind turbulence. 
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such a discrimination could be done. In fact, it is not even clear that interstellar 
turbulence can be described as an ensemble of Alfven waves. 



6.2 Do We Understand the Flat RM Structure Functions? 

Rotation Measure structure functions should have a logarithmic slope of 5/3 if 
both density and magnetic field fluctuations have a Kolmogorov spectrum, and 
an observationally indistinguishable value of 3/2 if a Kraichnan spectrum applies. 
This results holds if the two lines of sight are separated by a distance which is 
in the inertial subrange of the turbulence. Such a slope is rarely measured; a 
number of independent investigations have found that the RM structure functions 
on angular lags of several tenths of a degree to several degrees have logarithmic 
slopes of ~ 1/2, or even flatter. 

The interpretation of this result has been that the angular lags probed corre- 
spond to large scales in the interstellar medium, of order the outer scale or l arger 
(JMinter fc Spangler ll996tlHaverkorn et al ll20od . l2008l) . lMi"nter fe Spangler I (f~19961 



suggest that the aforementioned data are probing 2D turbulence that exists in 
sheets, and that the fully 3D component of the turbulence is on spatial scales less 
than the thickness of these sheets, with corresponding angular scales of order a 
few tenths of a degree or less. These analyses are, in fact, the basis of the claim 
that the outer scale is of order 1-5 parsecs in extent. 

Nonetheless, it would be comforting to actually measure, in a clear and un- 
ambiguous fashion, the transition from a ~ 1/2 logarithmic slope on angular lags 
> 1° to ~ 5/3 on angular lags < 0.2°. This would securely establish the value 
of the outer scale in the WIM turbulence. Without such measurements, we will 
continue to be tormented by the specter of on ISM in which we are measuring the 
inertial range of the turbulence (assuming this to be a meaningful concept), and 
that we lack an explanation for its logarithmic slope. 



6.3 What is the Significance of the "Pulsar Arcs" 

One of the most intriguing developments of the last decade in the study of in- 
terstellar turbulen ce has been the discovery of "pulsar arcs" ( Walker et al Il2004t 



ICordes et al 1 .2006) . The interpretation of this phenomenon seems to require that 
the turbulence responsible for interstellar turbulence is confined to one, or at most, 
a few thin sheets. A recent ov erview of the ob servational properties of the arcs and 
their interpretation is given in lRickett I ( 20111) . Is interstellar turbulence really con- 



fined to thin sheets? If so, what is the physical mechanism for the formation of 
them? If not, how are the Pulsar Arcs formed? 



6.4 What Generates Interstellar Turbulence? 

Plasma turbulence, as revealed by the density fluctuations responsible for radio 
scintillations, appears to be very widely distributed in the interstellar medium. 
The question then arises as to the mechanism responsible for its generation. 
A general consensus holds that the free energy source is expanding supernova 
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remnants, stellar superbubbles, and expanding HII regions ( Norman and Ferrara I 
Il996t iMacLow and Klessen I l2004t Iffill et al 1120121) . The small scale fluctuations 
that we detect in radio scintillations might be generated by baroclinic effectQ at 
the expanding interfaces between supernova remnants, stellar bubbles, and the 
ISM. However, there may be a problem with the distribution or diffusion of turbu- 
lence throughout the ISM. Supernova remnants and stellar bubbles occupy a very 
small fraction of the ISM, and turbulent damping limits the extent to which tur- 
bulence can propagate from the gener ation site to locations throughout the WIM 



(jSpangler 1120071: ISpangler et al Il201lh . 



6.5 Removal of Fast Magnetosonic Waves from Interstellar Turbulence 

Fast Magnetosonic waves are one of the three MHD wave modes, so one would 
expect them to comprise part of interstellar turbulence. However, it has been ar- 
gued that they can only constitute a negligibly small fraction of the energy in 
interstellar turbulence. The argument is based on the rapid damping of Fast Mode 
waves on thermal ions for conditions appropriate to the Warm Ionized Medium. 
If a sizeable fraction of the energy in interstellar turbulence is in the form of Fast 
Mode waves, then the large power input t o the interstella r med ium would ex- 
ceed the cooling capacity of the WIM gas ( Spangler1ll99ll I2005J) . Interestingly 



enoug h, they also seem to be abs ent from the solar corona and the solar wind at 
1 a.u.. lHarmon and Coles I (2005J) make a convincing argument that a substantial 



contribution of Fast Mode waves to the coronal tu rbulence budget is incompati- 
ble with spaced-receiver propagation measurements. iKlein et al I ( 20121) also argue 



that Fast Modes waves, or fluctuations possessing Fast Mode properties, consti- 
tute an insignifica nt portion of solar wind turbulence at 1 a.u.. The analysis of 
IKlein et al I (2012) is based on calculations of simulated turbulence, consisting of 



a superposition of Fast Mode waves, Slow Mode waves, and Alfven waves. These 
simulated realizations of turbulence are compared with actual measurements of 
solar wind turbul ence, especially the density-magnetic field correlation function. 
IKlein et al I (2012J) find that the realizations that resemble the true, observed tur- 



bulence are those with an insignificant fraction of Fast Mode waves. It should be 
mentioned before leaving this topic that the absence of Fast Mode waves in helio- 
spheric plasmas is a characteristic of plasmas far from shocks or other sources of 
unstable particle distributions. Shocks produce ion streaming instabilities which, 
in turn, generate beautiful, large amplitude Fast Magnetosonic waves, the best 
known examples o f which are the waves upstream of the Earth's bow shock (see 
iHoppe et all 1 198 ll for an entry point to a large literature). However, it seems to 
be the case that these Fast Mode waves are confined to relatively thin layers that 
bound strong shocks in the solar wind. To conclude this subsection, whether such 
a minor role for the Fast Mode in astrophysical turbulence is due to enhanced 
damping, or the turbulence generation mechanisms remains to be determined by 
future research. 



4 Baroclinic effects involve the generation of fluid vorticity by misaligned gradients of pres- 
sure and density. 
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6.6 Do We Understand the Lack of a Spectral Break at the Ion-Neutral 
Collisional Scale? 

As discussed in Section 3.1, there is observational evidence for a power law spec- 
trum of interstellar turbulence from scales of at least a few parsecs, down to scales 
as small as 100 km. The existence of a power law spectrum seems to indicate, 
on general grounds of dimensional analysis, that there are no fundamental scales 
between the outer scale on which stirring is done, and the inner scale where dissi- 
pation occurs. This assumption is in stark contrast to the situation for the WIM 
phase of the interstellar medium, in which a fundamentally-defined scale, the col- 
lisional scale l c = ^?S where Va is the Alfven speed and v in is the ion-neutral 
collisional scale, lies between the inner and outer scale s (l c — 10 15 — 10 16 cm). This 
conundrum was discussed by lArmstrong et al I (J1995I) . 

ICho and Lazarian I (|2003l this paper also co n tains references to earlier re- 
sults by these authors) and lOishi and MacLow I (J20061) present numerical sim- 
ulations showing that irregularities exist on scales smaller than the collisional 
scale, and conclude that ion-neutral collisional effects do not truncate the tur- 
bulent cascade. A different conclusion apears to be reached bv IShaikh and Zank I 
(|2008l ). IShaikh and Zank I (|2008l) find that while the turbulent cascade continues 
at wavenumbers larger than that corresponding to the ion-neutral collisional scale, 
the fundamental physics of the nonlinear interaction is modified. These authors 
claim that the magnetic field and velocity power spectra of the ionized fluid are sig- 
nificantly steepened in comparison with the same spectra of a fully- ionized plasma. 
This result would seem to be discordant with the observed Kolmogorov spectrum 
of density fluctuations down to much smaller spatial scales. This topic deserves 
further attention by the interested community. 



6.7 The Outer Scale of Interstellar Turbulence and Cosmic Ray Propagation 

In Section 3.6 above, we noted that several independent investigations find that 
the outer scale of turbulence in the WIM must be of orde r a few parsecs. There is 
an associated curiosity that was raised in lSpangler I (J200ll) , but does not appear to 
have been discussed since. If there is a break in the turbulence corresponding to an 
outer scale of 4 parsecs, then there should be an associated change in the transport 
properties of cosmic rays which are resonant with such irregularities, i.e. those with 
energies of 10 15 — 10 16 eV. In fact, the famous "knee" in the cosmic ray spectrum 
occurs here, but other mechanisms are normally invoked for its existence. Should 
the change in turbulence properties at scales that resonate with such cosmic rays 
be added to the mechanisms considered? In considering this matter, it should be 
recognized that the outer scale of turbulence in the Galactic halo is probably large, 
of order 100 pc (see Section 3.6 above). The larger fluctuations in the Galactic halo 
could resonate with higher energy cosmic rays than the smaller fluctuations in the 
WIM of the spiral arms. The fluctuations in the halo may dominate the Galactic 
transport of cosmic rays. 
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6.8 Can Observations Provide a Connection to Theories of Kinetic Processes in 
Turbulence? 

One of the most intriguing recent developments in the study of plasma turbu- 
lence has been elucidation of the role of kinetic processes in turbulence (i.e. those 
described by the Vlasov equation rather than MHD), and the observational sup- 
port for these ideas in spacecraft measurements of solar wind turbulence ( see 
lAlexandrova et al Il2009t and Section T3. 51 above') . Can we provide similar evidence 
of kinetic processes in the interstellar medium? Do the same kinetic processes 
which appear crucial in the solar wind, such as Landau damping of kinetic Alfven 
waves, play an important role in the interstellar medium? In kinetic damping pro- 
cesses, energy will flow to either ions or electrons, depending on which is resonant 
with the fluctuations being damped. However, observations of the WIM plasma, 
the best diagnosed astrophysical plasma, show temperature equilibration between 
electrons and different ion species (JHaffner et al 1 12009). A similar situation occurs 
in the clouds of the Very Local Interstellar Medium, where the same tempera- 
ture characterizes neutral atoms as well as several ion s with different masses (and 
therefore cyclotron frequencies) ( Spangler et al 11201 ll) . 



6.9 Why is the Spectrum of Plasma Turbulence the Same in HII Regions, the 
Warm Ionized Medium, and the Solar Wind? 

The I Armstrong et al I (|l995l) result of a Kolmogorov density fluctuation spectrum 
over several decades pertains to the WIM component of the interstellar medium, 
having been established from observations of relatively nearby pulsars and ex- 
tragalactic radio sources whose lines of sight are at high galactic latitudes. Radio 
wave propagation measur ements made on heavily-scattered lines of sight t hat pass 
through HII regions (e.g. ISpangler and Cordes Ill988t iMolnar et allll995l) are also 
consistent with a Kolmogorov density spectrum. Finally, the plasma of the so- 
lar corona and inner solar wind also has a Kolmogorov spectrum, particularly in 
the slow solar wind. The spectra of magnetic field and flow velocity also possess 
inertial subranges with spectra that are close to Kolmogorov, although recent re- 
search points to a slig ht difference in the power law indices of the magnetic field 
and velocity spectra ( Boldvrev et al 1120111) . This result is, perhaps, somewhat un- 



expected since the mechanisms for generation of the turbulence at the outer scale 
are presumably quite different in these different media, as might be the plasma /3 
that determines the dissipation mechanisms. Is the similarity of the turbulence in 
these quite different media a consequence of the universality turbulence? 



7 Summary and Conclusions 

Radio propagation observations yield a surprising amount of quantitative infor- 
mation about the plasma state of the interstellar medium, particularly turbulence 
in the WIM phase of the ISM and HII regions. The measurements emphasized 
in this paper have been Faraday rotation of linearly polarized signals that have 
propagated through the ISM during the passage from extragalactic radio sources 
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to the Earth, frequency-dependent polarization characteristics of the Galactic syn- 
chrotron radiation, and scintillations of Galactic and extragalactic radio sources 
due to small scale density fluctuations in the ISM. From this information, we can 
deduce the amplitude and spectral properties of interstellar turbulence. We have 
information on the outer and inner scales of this turbulence; the corresponding in- 
ertial subrange extends over roughly 10 decades. In some respects, interstellar tur- 
bulence resembles the extensively studied turbulence in the solar wind. However, 
there appear to be significant differences as well. In spite of impressive progress in 
this field, there are several (at least) poorly-understood aspects of interstellar tur- 
bulence that warrant the term "mysteries" . Several of these aspects are discussed 
in Section 6 above. Some of these could be addressed in a significant way with new 
observations on new or substantially upgraded radio telescopes such as the VLBA 
and LOFAR. 
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